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Abiotic  and  biotic  reductive  dechlorination  with  chlorinated 
ethene  dense  non-aqueous-phase  liquid  (DNAPL)  source  zones 
can  lead  to  significant  fluxes  of  complete  and  incomplete 
transformation  products.  Accurate  assessment  of  in  situ  rates 
of  transformation  and  the  potential  for  product  sequestration 
requires  knowledge  of  the  distribution  of  these  products  among 
the  solid,  aqueous,  and  organic  liquid  phases  present  within 
the  source  zone.  Here  we  consider  the  fluid— fluid  partitioning 
of  two  of  the  most  common  incomplete  transformation 
products,  c/s-1,2-dichloroethene  (c/s-DCE)  and  vinyl  chloride 
(VC).  The  distributions  of  c/s-DCE  and  VC  between  the  aqueous 
phase  and  tetrachloroethene  (PCE)  and  trichloroethene  (TCE) 
DNAPLs,  respectively,  were  quantified  at  22  °C  for  the 
environmentally  relevant,  dilute  range.  The  results  suggest 
that  partition  coefficients  (concentration  basis)  for  VC  and  c/s- 
DCE  are  70  ±  1  Laq/LjcE  dnapl  and  105  ±  1  Laq/LpcE  dnapl, 
respectively.  VC  partitioning  data  (in  the  dilute  region)  were 
reasonably  approximated  using  the  Raoult's  law  analogy  for 
liquid— liquid  equilibrium.  In  contrast,  data  for  the  partitioning  of 
c/s-DCE  were  well  described  only  when  well-parametrized 
models  for  the  excess  Gibbs  free  energy  were  employed.  In 
addition,  available  vapor— liquid  and  liquid— liquid  data 
were  employed  with  our  measurements  to  assess  the  temperate  re 
dependence  of  the  c/s-DCE  and  VC  partition  coefficients. 
Overall,  the  results  suggest  that  there  is  a  strong  thermodynamic 
driving  force  for  the  reversible  sequestration  of  c/s-DC  and 
VC  within  DNAPL  source  zones.  Implications  of  this  partitioning 
include  retardation  during  transport  and  underestimation  of 
the  transformation  rates  observed  through  analysis  of  aqueous- 
phase  samples. 

Introduction 

The  potential  of  reductive  dechlorination  processes  to 
transform  chlorinated  contaminants  to  henign  (or  less  toxic) 
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compounds  has  attracted  considerable  interest  within  the 
subsurface  remediation  community.  Source  zones  containing 
dense  non-aqueous-phase  liquids  (DNAPLs)  are  of  particular 
interest  il—3).  Engineering  enhancements  made  within 
DNAPL  source  zones,  such  as  the  supply  of  reactive  metal 
particles  or  the  stimulation  of  indigenous  microbial  popula¬ 
tions,  may  provide  an  effective  remediation  strategy  (e.g., 
refs  4  and  5).  Abiotic  and  biotic  reductive  dechlorination  of 
chlorinated  ethenes  can  proceed  via  hydrogenolysis,  a-e- 
limination,  /3-elimination,  hydrogenation,  and  coupling 
pathways  (6—8).  While  appreciable  quantities  of  coupling 
products  were  reported  in  some  metal-limited  systems  (e.g., 
refs  9  and  10),  heterogeneous  and  catalyzed  reactions 
involving  metals  generally  result  in  relatively  rapid,  complete 
dechlorination  to  ethene  and  ethane  via  a  combination  of 
/i- elimination,  hydrogenation,  and  hydrogenolysis  reactions 
(6).  Metabolic  reductive  dechlorination,  however,  proceeds 
via  a  slower,  sequential  hydrogenolysis  pathway  that  can 
result  in  the  accumulation  of  potentially  more  toxic  chlo¬ 
rinated  degradation  products  (S). 

Transformation  of  a  chlorinated  ethene  to  ethene  via 
metabolic  reductive  dechlorination  requires  a  consortium 
of  organisms,  with  complete  transformation  to  ethene  linked 
to  the  presence  and  activity  of  Dehalococcoides  spp.  (5,  7). 
Differences  in  electron  donor  availability  and  utilization, 
chlorinated  ethene  toxicity  limits,  and  pH  sensitivity  among 
the  multiple  species  required  for  the  transformation  of 
tetrachloroethene  (PCE)  and  trichloroethene  (TCE)  to  ethene 
are  just  some  of  the  challenges  to  the  effective  engineering 
and  management  of  this  dechlorination  strategy  (3, 5) .  When 
the  overall  dechlorination  process  stalls  (or  fails  to  completely 
dechlorinate  PCE  and  TCE),  c/s-dichloroethene  (c/5-DCE)  and 
vinyl  chloride  (VC)  accumulate  at  levels  well  above  regulated 
concentrations  (3, 11).  Given  that  the  production  of  c/s-DCE 
and  VC  corresponds  to  the  same  molar  extent  as  the  amount 
of  degradation  of  PCE  (or  TCE)  and  the  presence  of  DNAPLs 
(even  at  low  saturation)  represents  considerable  mass, 
degradation  products  can  approach  mUlimolar  concentra¬ 
tions  (11). 

While  accumulation  of  c/s-DCE  and  VC  during  DNAPL 
source  zone  treatment  is  strongly  influenced  by  microbial 
transformation  rates,  it  is  important  to  recognize  that  the 
production  of  c/s-DCE  and  VC  is  occurring  vnthin  multiliquid 
environments.  Thus,  remediation  design  requires  improved 
knowledge  of  the  interplay  between  microbial  transformation 
and  degradation  product  partitioning  (12,  13).  In  fact, 
degradation  product  partitioning  may  help  explain  mass 
balance  and  reactive  transport  anomalies  reported  in  labo¬ 
ratory-scale,  metabolic  reductive  dechlorination  studies  (e.g., 
refs  13—19).  Thus,  the  overall  objective  of  this  work  is  to 
elucidate  the  equilibrium  distribution  of  VC  and  c/s-DCE 
between  the  aqueous  phase  and  chloroethene  DNAPLs  across 
a  range  of  temperatures  relevant  to  source  zone  remediation. 
The  results  reported  here  are  used  to  illustrate  the  implica¬ 
tions  of  degradation  product  partitioning  within  DNAPL 
source  zones. 

Materials  and  Methods 

Chemicals.  Vinyl  chloride  (gas,  99.5-1-%  purity),  c/s- 1,2- 
dichloroethene  (97%  purity),  and  tetrachloroethene  (99% 
purity)  were  purchased  from  Sigma-Aldrich  and  used  as 
received.  Trichloroethene  (99-l-%  purity)  was  purchased  from 
Acros  Organics  and  used  as  received.  Dimethyl  sulfoxide 
(DMSO)  (99.9%  purity)  was  purchased  from  Fisher  Scientific. 
DMSO  was  used  to  minimize  volatilization  and  to  ensure 
miscibility  during  sample  preparation  and  analysis.  Water 
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was  purified  using  a  Millipore  Milli-Q  Gradient  AlO  system 
to  a  resistivity  of  >  18.2  MQ  •  cm  and  total  organic  carbon  of 
<15  ppb. 

Batch  Methods.  Batch  reactors  containing  Milli-Q  water, 
TCE,  and  VC  and  those  containing  Milli-Q  water,  PCE,  and 
ci.s-DCE  were  prepared  in  triplicate  using  Kimble  30  mL  glass 
centrifuge  tubes  sealed  with  poly  (tetrafluoroethene)  (PTFE)  / 
silicone  septa  and  open-top  screw  caps.  Because  it  was  not 
practical  to  directly  add  the  necessary  mass  of  VC  to  the 
reactors  in  gaseous  form,  a  sparging  procedure  was  developed 
to  create  solutions  of  VC  in  TCE  (Figure  SI-1,  Supporting 
Information) .  Given  the  carcinogenicity  of  VC,  a  management 
strategy  was  developed  and  implemented  to  ensure  the  safety 
of  all  personnel  (see  the  Safety  and  VC  Management  section 
in  the  Supporting  Information) .  Prior  to  using  the  VC-sparged 
TCE,  the  VC  concentration  was  quan tided  in  the  TCE  solution 
using  the  analytical  methods  described  below. 

Reactors,  containing  approximately  20  mL  of  MOli-Q  water 
and  10  mL  of  non-aqueous-phase  liquid  (NAPL)  comprising 
various  amounts  of  TCE  and  VC  or  PCE  and  ds-DCE,  were 
equilibrated  for  a  minimum  of  72  h  on  LabQuake  oscillating 
shaker  trays  at  22  ±  2  °C  (VC  studies)  or  22  ±  0.1  °C  (cd-DCE 
studies).  Previous  studies  have  determined  this  duration 
sufficient  to  reach  equilibrium  ( 1 6,  20) .  The  difference  in  the 
level  of  temperature  control  employed  in  these  studies 
resulted  from  safety  precautions  taken  to  ensure  VC  was 
contained  within  a  ventilating  fume  hood.  Reactors  contain¬ 
ing  c!5-DCE  were  allowed  to  mbc  within  a  constant-temper¬ 
ature  room  containing  a  recirculating  fume  hood.  Following 
the  period  of  equilibration,  the  tubes  were  centrifuged  at 
1500  rpm  for  10  min  to  ensure  separation  of  the  aqueous 
phase  and  organic  phases  (20).  Aqueous-phase  and  NAPL 
samples  were  collected  in  triplicate  for  subsequent  analysis. 
Subsequent  to  sampling  each  reactor  for  composition, 
additional  samples  of  both  phases  were  collected  for 
quantification  of  density.  Densities  of  the  aqueous  and 
organic  phases  were  determined  in  triplicate  using  2  mL 
glass  pycnometers  calibrated  daily  using  Milli-Q  water  at  22 
±  2  °C. 

Solubility  Measurements.  The  aqueous  solubility  of  ds- 
DCE  was  determined  as  described  above  with  reactors 
containing  only  water  and  excess,  neat  ds-DCE.  The  aqueous 
solubility  of  VC  could  not  be  determined  in  this  manner. 
Thus,  VC  gas  was  sparged  through  MUli-Q  water  at  22  ±  2 
°C  for  no  less  than  30  min  using  the  previously  described 
apparatus  (Figure  SI-1,  Supporting  Information).  Samples 
of  the  aqueous  phase  were  then  collected  in  quadruplicate. 
The  sampling  was  repeated  twice  more  after  two  10  min 
periods  of  sparging  to  produce  quadruplicate  samples  at  three 
sparging  durations  (or  12  total  samples). 

Analytical  Methods.  Triplicate  samples  of  the  aqueous 
phase  (20—70  uL)  and  NAPL  (20—50  uL)  were  collected  from 
the  reactors  containing  VC  using  a  gastight  syringe  and 
transferred  to  sealed  headspace  vials  for  subsequent  analysis 
as  described  below.  NAPL  samples  were  diluted  into  1.8  mL 
of  DMSO  and  then  subsampled  (10-50  uL)  in  triplicate  to 
headspace  vials  for  quantification.  Dilutions  were  varied  to 
be  consistent  with  the  linear  response  range  of  the  analytical 
method.  All  VC  samples  and  standards  contained  ~10  mL 
of22%  (w/w)  DMSO  in  Milli-Q  water.  For  reactors  containing 
cis-DCE,  triplicate  samples  were  collected  from  the  aqueous 
phase  (~5  mL)  and  NAPL  (~20 /<L).  Aqueous-phase  samples 
were  diluted  with  22  wt  %  DMSO  in  water  to  a  total  mass  of 
10  g.  NAPL  samples  were  diluted  into  8  mL  of  DMSO,  which 
was  then  subsampled  and  diluted  with  water  (10  x  dilution 
to  a  total  volume  of  10  mL)  to  produce  a  sample  for  analysis 
that  contained  ~11  wt  %  DMSO. 

VC  and  cis-DCE  concentrations  were  measured  using  a 
Hewlett-Packard  5890  series  II  gas  chromatograph  equipped 
with  an  HP  Plot  Q  column  (30  m  x  0.530  mm),  a  flame 


ionization  detector,  and  a  Perkin-Elmer  Turbo  Matrix  40  Trap 
headspace  sampler.  Headspace  vials  were  heated  to  95  °C, 
shaken  for  24  min,  and  pressurized  to  25  psi  (via  N2  addition) 
by  the  headspace  sampler  prior  to  transfer  of  3  mL  of  gas  to 
the  gas  chromatograph  for  analysis.  The  instrument  was 
calibrated  prior  to  each  day  of  use  with  a  six-point  calibration 
curve. 

Thermodynamic  Modeling.  Partitioning  of  contaminants 
between  aqueous  and  nonaqueous  phases  in  the  environ¬ 
ment  is  frequently  assumed  to  follow  an  analogy  to  Raoult’s 
law,  where  the  aqueous  solubility  of  the  contaminant  is 
multiplied  by  the  mole  fraction  of  the  contaminant  in  the 
nonaqueous  phase  to  produce  an  effective  solubility  (21): 


_  ^q.soly^APL 


(1) 


where  xf'  and  are  the  mole  fractions  of  component  i 
in  the  aqueous  phase  and  NAPL,  respectively,  and  is 
the  mole  fraction  representing  the  aqueous  solubility  of 
component  i.  Liquid— liquid  equilibrium,  however,  is  rigor¬ 
ously  described  by  equating  the  fugacities  of  each  liquid  phase 
(22): 


where  yf'i  and  y  are  the  activity  coefficients  of  component 
i  in  the  aqueous  phase  and  NAPL,  respectively,  0“  is  the 
fugacity  coefficient  of  component  i  at  the  vapor  pressure  P? 
(the  reference  state  for  component  i),  is  the  liquid-phase, 
partial  molar  volume  of  component  i,  R  is  the  universal  gas 
constant,  and  Tis  the  temperature  (K) .  It  should  be  recognized 
that  the  activity  coefficient  in  each  phase  is  a  function  of  x, 
and  that  the  Poynting  correction  (i.e.,  the  exponential 
function)  corrects  for  the  weak  pressure  dependence  of  the 
activity  coefficient. 

If  the  aqueous-phase  mole  fraction  of  component  i  falls 
within  the  dilute  solution  approximation,  then  yf’i  can  be 
replaced  by  the  infinite  dilution  activity  coefficient,  A 
similar  substitution  can  be  made  for  yf-^^  if  the 

values  of  x?"”^  permit  a  dilute  solution  approximation.  Under 
these  conditions,  the  liquid-phase  molar  volume  of  the  pure 
component,  ift,  may  be  used  to  approximate  (22).  Thus, 
these  conditions  produce  a  ratio  of  mole  fractions  that  is 
constant  and  equal  to  the  partition  coefficient  (on  a  mole 
fraction  basis)  for  the  dilute  region,  teq  3): 

^APL  ^  ^«,aq 

yAq  ~  .  ~.NAPL 
-'■i  /  i 

In  systems  where  the  dilute  solution  approximation  holds 
for  the  aqueous  phase  but  is  inappropriate  (or  questionable) 
for  the  NAPL,  Henry’s  law  may  be  employed  to  describe  the 
aqueous-phase  composition  in  equilibrium  with  a  hypo¬ 
thetical  real  gas  (22): 


-  ■Pref)'! 

<PiPi  =  exp^^ - — - j  (4) 

In  eq  4,  4>i  is  the  gas-phase  fugacity  coefficient  of 
component  i,  pi  is  the  partial  pressure  of  component  i  in 
the  gas  phase,  and  Pref  is  the  reference  pressure  at  which 
the  Henry  law  coefficient,  is  dehned  (often  taken  to  be 
Pi).  The  weak  pressure  dependence  of  iT,  is  accounted  for 
using  the  Poynting  correction,  again  assuming  that  p  is 
an  adequate  approximation  for  at  infinite  dilution.  For 
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systems  of  m  components  the  Henry  law  coefficient  for 
component  i  in  the  mixture,  is  dehned  as  (22) 


m  m—l  m 

In  ^  Xj  in  H^j  -  X  X  t5) 

j=i  y=i  k>j 

where  fl,i-  are  interaction  parameters  between  the  yth  and 
fcth  components.  Many  multicomponent  systems  of  en¬ 
vironmental  relevance  contain  dilute  solutes,  so  H,,niix  can 
he  well  approximated  as  the  Henry’s  coefficient  for  the 
binary,  Hy  (e.g.,  water  and  volatile  contaminant),  where 
the  solute  and  solvent  are  components  land  j,  respectively. 

The  fugacity  coefficient  in  the  gas  phase  can  be  calculated 
using  the  second  virial  coefficient  and  a  Pitzer  expansion  as 
(22) 


ln0,=  ^2Xy&(4+.,.4)- 


with 


4- 


c,ij  j 

(6a) 

im 

(6b) 

iW) 

(6c) 

where  y,  and  y,  are  the  mole  fractions  of  components  i  and 
j  in  the  gas  phase,  T  is  the  temperature,  and  Tc  y  and  Pc.y  are 
the  critical  temperature  and  pressure,  respectively,  of  the 
mixture  formed  by  components  i  and  ].  Tc.y  and  Pcjj  are 
dehned  using  standard  mixing  rules  (22). 

Use  of  Henry’s  law  (eq  4)  to  describe  the  aqueous-phase 
fugacity  in  eq  2  produces  the  following  equation,  which 
describes  the  partitioning  of  a  solute  that  is  assumed  to  be 
dilute  within  the  aqueous  phase: 


While  is  often  a  function  of  the  composition  (20),  it  can 
become  constant  when  is  not  a  function  (or  is  a  weak 
function)  of  the  NAPL  composition,  usually  under  an  ideal 
solution  assumption  (yNAP^  =  i).  if  the  NAPL  is  assumed  to 
be  an  ideal  solution  (i.e.,  =  1),  eq  7  can  be  simplihed 

to  produce  the  Raoult  law  analogy  (RLA)  (eq  8)  for  partihoning 
of  component  (between  a  dilute  aqueous  phase  and  an  ideal 
NAPL. 


Comparison  of  eqs  1  and  8  dehnes  inverse  of 

It  should  be  noted,  however,  that  the  reference  state 
for  component  i  is  treated  as  a  real  gas;  thus,  Raoult’s  law 
is  never  strichy  applied.  The  RLA  simply  represents  an  analogy 
between  eq  1  and  that  of  Raoult’s  law.  Alternatively,  the 
partition  coefficient  can  be  (i)  quantified  by  regressing 
liquid— liquid  equilibrium  data  at  a  specific  temperature  or 
(ii)  estimated  for  any  set  of  conditions  (including  the  limiting 
conditions  that  produce  and  using  a  group  con¬ 


tribution  method  (e.g.,  UNIFAC)  to  simulate  liquid  activity 
coefficients  for  use  in  an  isothermal  flash  calculation. 

Use  of  the  RLA  to  calculate  the  partition  coefficient  implies 
that  the  Henry  law  coefficient  is  known  at  all  temperatures 
of  interest.  The  temperature  dependence  of  Henry’s  law 
coefficients  is  frequently  modeled  using  the  Van’t  Hoff 
equation,  which  performs  well  over  a  small  temperature  range 
where  the  partial  molar  enthalpy  of  volatilization  can  be 
assumed  to  be  constant.  This  assumption,  however,  fails  for 
most  solutes  over  an  expanded  range  of  temperatures, 
resulting  in  the  use  of  empirical  descriptions  for  the  tem¬ 
perature  dependence  of  Henry’s  law  coefficients  (22,  23). 
Recently,  an  empirical  model  (eq  9)  was  employed  to  describe 
the  Henry  law  coefficients  for  organic  solutes  in  water  by 
assuming  that  the  partial  molar  enthalpy  of  volatilization  is 
a  linear  function  of  the  temperature  (23). 

H,  =  exp(A  -  I  -f  Cln  t]  (9) 

A,  B,  and  C  in  eq  9  are  htting  parameters.  Regression  of  eq 
9  to  existing  H,  data  was  accomplished  using  a  weighted 
least-squares  approach  (24).  A  10%  relative  standard  error 
was  assumed  for  all  data  sets  not  reporting  uncertainty. 

Results  and  Discussion 

Use  of  the  RLA  requires  reliable  values  of  the  aqueous 
solubility.  Thus,  existing  data  were  reviewed  and  reanalyzed 
to  develop  new  correlations  for  the  temperature  dependence 
of  the  Henry  law  coefficient  (and  hence  solubility)  describing 
the  partitioning  of  VC  and  cw-DCE  between  aqueous  and 
gas  phases.  Data  are  available  in  the  form  of  vapor— liquid 
equilibria  or  liquid— liquid  equilibria  {13, 25—34),  EPICS  data 
(35,  36),  and  infinite  dilution  activity  coefficients  (37,  38). 
Many  of  these  vapor— liquid  equilibria  and  EPICS  studies, 
however,  overlook  nonidealities  in  the  gas  phase  that  are 
known  to  exist  for  VC  and  to  a  lesser  extent  cfs-DCE  (25)  and 
were  therefore  reanalyzed  using  eq  4.  The  EPICS  equation 
for  a  nonideal  gas  phase  (eq  SI-4,  Supporting  Information) 
was  developed  using  the  concentration  form  of  eq  4  and 
system  mass  balance  following  the  derivation  provided  by 
Gossett  (35). 

Regressions  were  accomplished  using  Hi  and  are  shown 
in  Figures  SI-2  and  SI-3  (Supporting  Information)  in  terms 
of  the  concentration  form  of  the  Henry  law  coefficient,  Hf, 
which  assumes  that  the  molar  density  of  the  solution  can  be 
approximated  by  the  molar  density  of  the  solvent  (i.e,  the 
solute  is  dilute  enough  to  have  a  negligible  influence  on  the 
physical  properties  of  the  solution).  This  assumption  was 
valid  for  VC  and  cIs-DCE  as  the  measured  densities  were 
within  0.2%  of  the  corresponding  density  of  water.  Values 
for  the  parameters  in  eq  9  were  found  to  be  A  =  214.8,  B  = 
11  688,  and  C=  -29.6  for  VC  {If  =  0.98)  and  A  =  182.0,  B 
=  11  158,  and  C  =  —24.4  for  ds-DCE  {If  =  0.86).  The  lower 
value  in  the  case  of  ds-DCE  results  from  the  large  variability 
in  the  measured  values  of  H,  over  a  limited  range  of 
temperature  (Figure  SI-3).  The  correlations  developed  here 
offer  the  improvements  of  greater  data  support  and  greater 
temperature  range  when  compared  to  previously  reported 
correlations  (35,  36).  Use  of  the  correlations  is  appropriate 
for  temperatures  of  0.2—85  °C  for  VC  and  10—40  °C  for  ds- 
DCE— an  important  advance  when  assessing  partitioning 
(gas— liquid,  solid-liquid,  or  liquid— liquid)  where  contami¬ 
nant  reduction  follows  thermal  remediation  (e.g.,  refs  39  and 
40). 

Aqueous  solubilities  of  VC  and  ds-DCE  (in  mole  fraction) 
calculated  using  eq  2  and  the  correlations  developed  for  the 
Henry  law  coefficients  (at  1  atm)  are  shown  as  a  function  of 
temperature  in  Figures  1  and  2.  Mole  fractions  were  converted 
to  concentration  by  assuming  that  the  presence  of  VC  or 
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- regressed  model  for  /-/,  at  1  atm 

- model  isobars  with  pressures  labeled  in  atm 

O  Hayduk  and  Laudie  (1 974)  adjusted  to  1  atm 
A  Carr  et  al.  (2000)  adjusted  to  1  atm 
V  Garrens  et  al,  (1967)  adjusted  to  1  atm 

□  DeLassus  and  Schmidt  (1981)  adjusted  to  1  atm 
O  Patel  et  al.  (1979)  adjusted  to  1  atm 

□  Gossett  (1987)  measured  at  1  atm 

O  Ashworth  et  al.  (1988)  measured  at  1  atm 
O  Berens  (1975)  adjusted  to  1  atm 

-  Wilhelm  (1977)  correlation  at  1  atm 

-  Berens  (1975)  correlation 

•  this  study  measured  at  1  atm 

o  Hayduk  and  Laudie  (1974)  measured 

*  Carr  et  al.  (2000)  measured 

»  Garrens  et  al.  (1967)  measured 
■  DeLassus  and  Schmidt  (1981)  measured 
o  Patel  et  al.  (1979)  measured 

♦  Berens  (1975)  measured 

FIGURE  1.  Solubility  of  vinyl  chloride  in  water.  Error  bars 
(standard  error)  are  included  for  data  at  1  atm  only.  Many  error 
bars  are  smaller  than  the  associated  symbol. 


regressed  model  for  H,. 

o 

McGovern  (1943) 

A 

Sato  and  Nakijima  (1979) 

□ 

Wright  et  al.  (1992) 

o 

Flick  (1998) 

□ 

Gossett  (1987) 

A 

Tseetal.  (1992) 

V 

Bissonette  et  al.  (1990) 

O 

Ashworth  et  al.  (1988) 

V 

Ervin  et  al.  (1980) 

— 

Hovarth  regression 

• 

this  study 

FIGURE  2.  Solubility  of  c/s-DCE  in  water.  Error  bars  represent 
the  standard  error. 

cw-DCE  had  a  negligible  influence  on  the  molar  volume  of 
the  aqueous  solution.  The  solubility  of  VC  was  available  at 
total  pressures  ranging  from  0.32  to  17.4  atm  (Figure  1) .  Each 
value  for  solubility  was  adjusted  using  the  Henry  law 
coefficient  at  1  atm  and  a  total  pressure  of  1  atm  for  a  binary 


O  VC  data 

- Rmodel  with  Raoult's  Law  analogy,/:!'-'^  =371 

-  regression  of  VC  data  =  348 


FIGURE  3.  Equilibrium  distribution  of  vinyl  chloride  between  an 
aqueous  phase  and  TCE  NAPL  at  22  ±  2  °C.  The  inset  provides 
better  visualization  of  xf  data  below  6  x  10  ^  Error  bars 
represent  the  standard  error  established  from  triplicate  sampling 
of  the  same  reactor. 

system  of  VC  and  water  (eq  2) .  In  the  case  of  VC,  the  pressure 
dependence  of  the  gas  solubility  is  shown  as  a  series  of  isobars 
in  Figure  1  and  contours  in  Figure  Sl-4  (Supporting  Infor¬ 
mation).  Also  shown  in  Figure  1  are  two  correlations  for  the 
aqueous  solubility  that  were  developed  in  the  1970s  by 
regressing  subsets  of  these  solubility  data.  The  new  correlation 
(in  the  form  of  a  Henry  law  coefficient)  better  captures  the 
large  data  set  and  offers  a  larger  range  of  applicable 
temperatures.  Comparison  of  the  measured  aqueous  solu¬ 
bility  of  VC  (49  ±  2  mM  at  22  ±  2  °C  and  1  atm)  to  that 
predicted  from  the  new  correlation  at  1  atm  and  22  °C  (42.4 
mM)  suggests  the  measurement  may  be  high.  Note  that  values 
of  Hi  calculated  using  these  measured  solubilities  are  shown 
in  Figures  SI-2  and  SI-3  (Supporting  Information),  but  the 
measurements  were  not  employed  in  the  regression.  It  is 
likely  that  the  discrepancy  between  the  measured  VC 
solubility  and  literature  values  resulted  from  a  pressure  in 
the  sparging  apparatus  that  was  slightly  greater  than  1  atm 
(headspace  pressures  were  not  monitored).  Also,  it  is 
important  to  recognize  that  the  assumption  of  1  atm  of  total 
pressure  also  means  that  there  is  no  NAPL  present  (note  the 
solubility  curve  at  the  vapor  pressure  in  Figure  1).  Any 
system  containing  NAPL  requires  use  of  the  aqueous  solubility 
corresponding  to  the  vapor  pressure  of  VC  (see  Table  SI-1, 
Supporting  Information,  for  Antoine  equation  parameters). 
The  measured  value  for  the  aqueous  solubility  of  ds-DCE 
(70  ±  4  mM  at  22  ±  2  °C)  appears  to  fall  within  the  scatter 
in  the  literature  values,  but  is  noted  to  be  higher  than  that 
predicted  by  the  correlation  at  22  °C  (60.5  mM). 

A  series  of  batch,  liquid— liquid  equilibrium  experiments 
were  employed  to  measure  the  distribution  of  VC  in  TCE 
NAPL  and  cfs-DCE  in  PCE  DNAPL.  Data  from  these  experi¬ 
ments  are  shown  in  Figures  3  and  4,  respectively,  in  terms 
of  mole  fractions.  The  insets  in  each  Hgure  show  the  lower 
concentration  region  of  each  plot.  Regressions  (linear  least- 
squares  weighted  by  the  relative  standard  error)  were  used 
to  determine  for  VC  (348  ±  4,  =  0.99)  and  cfs-DCE  (593 

±  5,  =  0.99).  These  regressions  assume  that  the  partition 

coefficient  was  not  a  strong  function  of  the  composition. 
The  regressed  values  for  the  partition  coefficients  for  VC 
between  an  aqueous  phase  and  TCE  DNAPL  and  for  ds-DCE 
between  an  aqueous  phase  and  PCE  DNAPL  can  be  ap- 
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FIGURE  4.  Equilibrium  distribution  of  cis-1,2-dicbloroetbene 
between  an  aqueous  phase  and  PCE  DNAPL  at  22  ±  0.1  °C.  Tbe 
inset  provides  better  visualization  of  xf  data  below  1.2  x  10  ^ 
Error  bars  represent  tbe  standard  error  established  from 
triplicate  sampling  of  the  same  reactor. 

proximated  on  a  concentration  basis  {K^‘)  by  70  ±  1  Laq/ 
Lice  DNAPL  ^nd  105  it  1  Laq/LpcEDNAPLi  respectively. 

Alternatively,  UNIFAC  [41,  42)  estimates  of  the  ratio  of 
infinite  dilution  activity  coefficients 

VC  and  cfs-DCE,  respectively  (not  plotted) .  It  is  hypothesized 
that  the  poor  performance  of  UNIFAC  when  estimating  the 
partitioning  of  VC  relates  to  the  use  of  compounds  with  lower 
vapor  pressures  («Fvc  =  3.65  atm  at  22  °C)  when  regressing 
the  interaction  parameters  between  groups  7  (water)  and  37 
(C=CC1)  (e.g.,  refs  42  and  43).  This  hypothesis  is  supported 
by  the  considerably  better  UNIFAC  estimate  of  the  partition 
coefficient  for  cfs-DCE,  which  is  within  10%  of  the  measured 
value.  In  addition,  this  hypothesis  suggests  that  UNIFAC  may 
be  employed  to  estimate  cfs-DCE  partitioning  in  TCE  NAPE. 

Application  of  the  RLA  produces  values  for  gf  37]^ 
and  917  for  VC  and  ds-DCE,  respectively  (Figures  3  and  4). 
The  RLA  estimate  of  the  VC  partition  coefficient  appears  to 
describe  the  data  well,  especially  if  data  that  fall  within  25 
X  10“^>x?'>>30  X  10“^  influenced  the  regression  to  produce 
a  lower  value  oUcj',  although  exclusion  of  these  data  from  the 
regression  produced  a  statistically  insignificant  change  to 
the  partition  coefficient  (p  =  0.34).  While  the  RLA  is  not 
capable  of  exacfly  reproducing  the  data,  it  does  provide  a 
reasonable  first  approximation  (±10%)  for  fhe  partitioning 
of  VC  between  the  aqueous  phase  and  NAPL.  Use  of  the  RLA 
(eq  8)  may  prove  advantageous  given  that  estimates  (on  a 
mole  fraction  basis)  are  independent  of  the  NAPL  composi¬ 
tion  (recall  that  the  RLA  presumes  the  NAPL  is  an  ideal 
solution).  Estimates  of  for  VC  over  the  range  of 

temperatures  (0.2—85  °C)  supporting  the  H,  model  are 
consolidated  using  the  following  equation: 


=  “2.04  X  10^  ±  15.3r- 

(2.43  X  10“^)7^  273.35  K  >  r>  358.15  K  (10) 

The  poor  performance  of  the  RLA  for  cis-DCE  may  be 
attributable  to  the  scatter  in  the  data  used  to  parametrize 
the  Hi  model  (Figure  Sl-3,  Supporting  Information).  Use  of 
the  Gossett  correlation  (which  produces  lower  values)  to 
estimate  H,  for  cfs-DCE  improves  the  RLA  estimate  (759) ,  but 
does  not  capture  the  data.  UNIFAC  performs  better  for  cis- 
DCE  than  for  VC,  providing  a  reasonable  first  approximation 
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FIGURE  5.  Reversible  sequestration  of  products  as  a  function  of 
NAPL  saturation:  fraction  of  mass  represented  by  tbe  aqueous- 
pbase  concentration  and  retardation  factor  due  to  partitioning. 


for  the  partition  coefficient.  UNIFAC  was  employed  with 
established  temperature  dependencies  of  fhe  mutual  solu¬ 
bilities  of  PCE  and  water  (44)  to  estimate  the  temperature 
dependence  of  the  partition  coefficient  for  cfs-DCE  in  PCE 
DNAPL  source  zones  (eq  11). 


^?;d.-DCE  =  6.659  X  10^  -  33.76r± 

(4.536  X  10“^) 273.15  K  >  r>  343.15  K  (11) 


These  new  correlations  (eqs  10  and  11)  along  with  the 
supporting  data  and  analysis  provide  additional  information 
necessary  to  better  understand  the  distribution  of  contami¬ 
nants  and  their  degradation  products  in  multiphase  sub¬ 
surface  environments  typical  of  those  found  at  chlorinated 
DNAPL  source  zones  for  a  variety  of  environmentally  relevant 
temperatures. 

Implications.  Partitioning  of  degradation  products  within 
DNAPL  source  zones  has  the  potential  to  reversibly  sequester 
product  mass.  Reversible  sequestration  may  have  two 
important  implications:  (i)  product  mass  located  within  the 
NAPL  will  likely  be  neglected  when  assessing  degradation 
kinetics  and  (ii)  the  transport  of  degradation  product 
downgradient  may  be  retarded  (due  to  partitioning).  These 
implications  are  illustrated  using  the  measured  partition 
coefficients  for  VC  and  cfs-DCE  to  calculate  the  fraction  of 
total  mass  that  is  located  within  the  aqueous  phase  (eq  SI-9, 
Supporting  Information)  and  the  retardation  factor  (eq  SI- 
10,  Supporting  Information)  for  a  range  of  NAPL  saturations 
assuming  local  equilibrium  (Figure  5).  If  partitioning  is  fast 
relative  to  degradation,  as  accumulation  of  ds-DCE  and  VC 
may  suggest,  then  laboratory  and  field  assessments  of 
degradation  kinetics  in  the  presence  of  NAPL  may  exclude 
large  fractions  of  mass  (e.g.,  aqueous-phase  concentrations 
account  for  <40%  of  the  total  product  mass  at  the  modest 
NAPL  saturation  of  2%). 

Several  authors  have  suggested  degradation  product 
partitioning  as  a  competing  process  that  influences  parenf 
and  producf  availability  and,  hence,  apparent  rates  of 
degradation  [13,  15—17,  45,  46).  Thus,  the  influence  of 
degradation  product  partitioning  on  the  metabolic  reduc¬ 
tive  dechlorination  process  occurring  within  a  PCE  DNAPL 
source  zone  was  assessed  using  illustrative  simulations. 
Here,  the  model  and  kinetic  parameters  described  by  Christ 
and  Abriola  (46)  are  employed  to  simulate  a  30  cm  column 
in  which  the  first  15  cm  contains  entrapped  PCE  DNAPL 
at  a  uniform  saturation  of  2%  (simulations  are  summarized 
in  the  Supporting  Information) .  Partitioning  was  observed 
to  produce  pseudosteady  effluent  concentrations  of  cis- 
DCE  and  PCE  that  are  approximately  40%  and  90%, 
respectively,  of  the  effluent  concentrations  produced  in 
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FIGURE  6.  Effluent  concentrations  (top)  and  pseudo-steady-state 
resident  aqueous-phase  concentrations  at  150  days  (bottom)  for 
simulations  conducted  with  (solid  line)  and  without  (dashed 
line)  degradation  product  partitioning.  The  column  is  30  cm 
long  with  a  PCE  DNAPL  source  zone  (s'^'  =  0.02)  located  in  the 
first  15  cm.  TCE  concentrations  (not  shown)  are  below  5  /iM. 

the  absence  of  partitioning  (Figure  6).  In  addition,  resident 
aqueous-phase  concentrations  (shown  at  150  days)  dem¬ 
onstrate  partitioning  attenuates  cis-DCE  concentrations 
within  the  15  cm  PCE  DNAPL  source  zone  (Figure  6) .  While 
the  interplay  between  partitioning  and  degradation  re¬ 
quires  further  study,  the  results  suggest  that  transformation 
rates  may  be  underestimated  and  remediation  systems 
may  be  underdesigned  when  degradation  product  parti¬ 
tioning  is  neglected. 

Recognition  that  reversible  sequestration  may  influence 
assessments  of  degradation  kinetics  requires  reassessment 
of  dosing  requirements  (both  the  total  amount  and  the  rate 
of  addition)  when  employing  soluble  substrate  to  produce 
the  required  electron  donor.  The  introduction  of  insoluble 
substrate  such  as  vegetable  oils  within  the  source  zone  (or 
immediately  downgradient)  may  only  accentuate  the  revers¬ 
ible  sequestration  [12,  15,  19).  Absorption  of  degradation 
products  (and  parent  contaminants)  by  insoluble  amend¬ 
ments,  while  temporally  decreasing  the  mass  discharge,  may, 
in  fact,  make  difficult  the  assessment  of  the  actual  kinetics 
and  extents  (contaminant  mass  transformed)  of  the  degra¬ 
dation  process.  Ongoing  efforts  are  aimed  at  elucidating 
partitioning  kinetics  in  NAPL  source  zones  so  that  the 
interplay  between  partitioning  and  degradation  can  be 
rigorously  assessed. 
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